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We propose a qualitatively new mechanism for quarkonium production, motivated by the global 
features of the experimental data and by the successes/failures of existing models. In QCD, heavy 
quarks are created in conjunction with a bremsstrahlung color field emitted by the colliding partons. 
We study the effects of perturbative gluon exchange between the quark pair and a comoving color 
field. Such scattering can flip the spin and color of the quarks to create a non- vanishing overlap with 
the wave function of physical quarkonium. Several observed features that are difficult to understand 
in current models find simple explanations. Transverse gluon exchange produces unpolarized J/ip's, 
the Xci and Xc2 states are produced at similar rates, and the anomalous dependence of the J/ip 
cross section on the nuclear target size can be qualitatively understood. 



I. INTRODUCTION AND SUMMARY 
A. Introduction 

Quarkonium production has turned out to be a chal- 
lenge as well as an inspiration for our understanding of 
hard QCD processes [Q-Q. In the case of standard in- 
clusive processes, the theoretical framework is uniquely 
defined by the QCD factorization theorem ||. This the- 
orem allows a physical cross section a to be expressed as 
a product of universal parton distribution and fragmen- 
tation functions multiplied by a subprocess cross section 
a, which is calculable in PQCD. The factorization theo- 
rem relies on a completeness sum over the final state and 
does not apply to the quarkonium cross section, which 
constitutes only a small fraction of the total heavy quark 
production cross section. 

While a theoretical description of quarkonium produc- 
tion is thus more model dependent, it can potentially re- 
veal more about the dynamics of hard processes than can 
be learned from, eg., the total heavy quark cross section. 
In particular, it is intuitively plausible that the quarko- 
nium cross section is sensitive to reinteractions with par- 
tons created along with the heavy quark pair. Thus 
quarkonium production can serve as a 'thermometer' of 
the environment, as has been recognized in the search for 
a quark-gluon plasma in heavy ion collisions In this 
paper we wish to explore the possibility that rescattering 
of the heavy quarks causes the puzzling anomalies seen 
in quarkonium hadroproduction. 



There is independent evidence that the environment 
in charm hadroproduction is rather 'hot'. In tt~ N — > 
DD + X the observed spread in relative azimuthal angle 
of the Z?-mesons requires an average intrinsic transverse 
momentum of the incoming partons (k^) ~ 1 GeV 2 JlO| . 
The 'leading particle' asymmetry between D~ and D + 
is larger than expected from PQCD, and persists for D- 
mesons produced with kj_ < 10 GeV 2 111. Both effects 
are weaker for photoproduced charm |l0| , |12j | . 

We shall study the effects of perturbative gluon ex- 
change between the heavy quark pair and a comoving 
Jl3| color field. The interaction is assumed to occur at 
an early stage, before the pair has expanded to the size 
of physical quarkonium. Hence only comovers which are 
created (via bremsstrahlung) in the hard process itself 
are relevant, whereas interactions with beam fragments 
at typical hadronic distances ~ 1 fm arc ignored. In this 
sense our approach differs from that of the 'Color Evap- 
oration Model' (CEM) @Jl|, which only considers late, 
non-perturbative interactions of the heavy quarks. On 
the other hand, similarly to the CEM our quark pairs 
are produced near threshold. Hence many of the phe- 
nomenological successes of the CEM concerning the de- 
pendence of quarkonium cross sections on various kine- 
matical variables are incorporated in our model. The 
'Color Octet Model' (COM) @, also considers late in- 
teractions, through an expansion in powers of the relative 
velocity v of the bound quarks as specified by NRQCD 
p7| . This expansion is general and should hold for any 
description of quarkonium production, including ours. 
The higher order (v/c) n terms need not, however, give 
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a dominant contribution to the cross section. To our 
knowledge, the COM assumption that the heavy quark 
pair is unaffected by earlier reinteractions with its envi- 
ronment has not been proven. 

Data on charmonium and bottomonium production is 
available for a wide variety of beams, targets and kinc- 
matical conditions. Comparisons with the COM and 
CEM approaches have met with some successes, but also 
with difficulties |l|-Q]. The data suggests a production 
dynamics which in some respects differs from the late 
and soft reinteraction scheme of the CEM and COM. In 
particular, 

(A) The heavy quark pair turns color singlet at an early 
stage, while the pair is still compact (ie., small com- 
pared to the size of the quarkonium wave function) . 

(B) In hadroproduction there is at least one secondary 
gluon exchange after the primary heavy quark pro- 
duction vertex. 

(C) The 'anomalies' of quarkonium production depend 
only weakly on the quark mass to, on the CM en- 
ergy and on the transverse momentum p± . 

We discuss the experimental basis for these features in 
Section || and then develop our QCD scenario in Sec- 
tion p| . This scenario applies to quarkonium production 
at both moderate and large transverse momentum. How- 
ever, we shall limit our discussion and the calculations 



/i 2 <C to 2 . We will investigate the effects of rescattering 
at this hardness scale fx. 



presented in Section IV to the total quarkonium cross 

In the rest of the 
Conclusions 



section, ie., moderate xf and p± ~ to. 
present Section we summarize our results, 
and an outlook are presented in Section [v|. 



B. Summary 

The basic Color Singlet Mechanism (CSM) JT|], which 
is known to grossly underestimate the J r /ip hadroproduc- 
tion cross section, is shown in Fig. la. The gluon emission 
takes place in a (proper) time r ~ 1/to, simultaneously 
with the heavy quark production process. This is com- 
patible with feature (A), but since there are no relevant 
later interactions the CSM does not agree with (B). The 
situation is qualitatively the same for loop corrections to 
the CSM (Fig. lb), since the space-time scale of the loop 
is 1/to. 

Prior to the heavy quark production vertex the collid- 
ing partons radiate gluons as part of the normal QCD 
structure function evolution. The space-time scale of 
this process is determined by the virtuality k 2 of the 
partons which couple to the heavy quark line. As is 
characteristic of evolution processes, the k 2 distribution 
a s (k 2 )dk 2 jk 2 is logarithmic between a lower cutoff de- 
termined by the (perturbative) factorization scale and an 
upper limit given by the heavy quark mass to. Thus the 
effective value of \k 2 \ is given by a perturbative scale which 
we denote by /A This scale grows with to 2 but satisfies 



^ 5>>3 ^ (la) 



JA|/ 





(lc) 

FIG. 1. Basic processes for J/ip hadroproduction in the 
CSM (figures (la) and (lb)) and in our model (figure (lc)). 



The approach presented in Section |I| is based on a 
perturbative reinteraction of momentum transfer t <~ fi 
between the heavy quarks and a classical color field T 
(Fig. lc). This field is assumed to originate from gluon 
bremsstrahlung in the gg — > QQ subprocess. The reinter- 
action can occur long after the heavy quarks are created 
provided the field T is comoving with the quark pair. 
Hence we shall assume T to be isotropic in the pair rest 
frame. The scale /i of the secondary interaction is smaller 
than the quark mass scale to of the CSM but larger than 
the bound state scale a s m of the COM0. The existence 
of the comoving color field T in hadroproduction is our 
main postulate, motivated by the data. There should be 
no corresponding field in the current fragmentation re- 
gion of photo- and leptoproduction, since photons do not 
radiate gluons (at lowest order). 

Our quarkonium hadroproduction amplitude is essen- 
tially given by the perturbative diagram of Fig. 2 (to- 
gether with the diagram corresponding to rescattering of 
the antiquark). Here $ is the (color octet) wave function 



*For the charmonium system, some of these scales are nu- 
merically similar, but should be distinguished for reasons of 
principle. The scales do differ for bottomonium. 
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k 2 -q 

FIG. 2. A perturbative interaction between the quark pair 
and a gluon from the color field V creates an overlap between 
the QQ wave function $ from the gg — > QQ process with 
the physical quarkonium wave function ^ . There is a second 
diagram where the gluon attaches to the antiquark. 

of the heavy quark pair produced in the fusion process 
99 ~ * QQ and ^ is the (color singlet) quarkonium wave 
function. 

As we shall see, there are two production mechanisms 
for spin triplet S'-wave quarkonia such as the J/V'- In 
the gg — > QQ subprocess the quark pair can either be 
produced in an S = L = state, followed by a spin-flip 
interaction with a transverse gluon from the color field T, 
or the pair can be produced with S = L = 1, followed by 
a spin-conserving interaction with a longitudinal gluon. 
The first contribution gives unpolarized quarkonia, since 
the quark pair is produced with total angular momentum 
J = and the color field V is isotropic. The second con- 
tribution turns out to give quarkonia with a transverse 
polarization. The striking experimental fact |l9|-|2l|] that 
the J /ip and t/j' are produced unpolarized (at moderate 
p±) thus implies that the former mechanism dominates, 
ie., that the gluons in T are transversely polarized. 

The P-wave quarkonium states xj are produced from 
quark pairs with S = 0, L — 1 followed by a spin-flip 
interaction with a transverse gluon from T. The calcu- 
lated cross section ratio c(Xi)/ (7 (X2) = 3/5, in agreement 
with the ratio measured in irN collisions for charmonium, 
0.6 ± 0.3 f^H]. We also find that due to the indirect 
XcJ ~ * J/"0+7 contributions the total J hp polarization is 
slightly longitudinal, A = —.14 (c/. Eq. (pq)). Taking into 
account also the CSM mechanism, which dominantly pro- 
duces Xc2's with J z = ±2 that decay into transversely po- 
larized J/ip's the expected ratio cr(xi) / <j(x2) < 3/5 
and A > —.14. This is compatible with data as long as 
the CSM does not dominate the rescattering mechanism 
(Fig. 2) for X2 production. 

The xj wave function VE" (Fig. 2) vanishes at the origin, 
which suppresses its overlap with small sized (~ 1/m) 
heavy quark pairs. Thus in the CSM the relative produc- 
tion rate of P- and S'-wave charmonium states is governed 
by the small ratio |i? x '(0)/2m| 2 /|i? v ,(0)| 2 ~ 0.01 ||. In 
our approach the initially compact quark pair expands, 
with velocity ~ fJ>/fn, before the gluon exchange in Fig. 2 
gives the quark pair the bound state quantum numbers. 
The wave function $ is thus an expanded version of the 



quark pair wave function created in the gg — > QQ sub- 
process. We model this by a scale factor p, which we fit to 
the measured a(xcj)/ '&dir{J /VO cross section ratio. We 
find that we need p ~ 3, suggesting significant expansion 
of the quark pair before reinteraction in XcJ production. 

Our approach can also be applied to quarkonium pro- 
duction at high ^ m, where the dominant production 
mechanism is gluon fragmentation |2q| . The fragmenting 
gluon is initially highly virtual and radiates gluons with 
hardness ranging from the factorization scale up to p±. 
The gluons of relatively small hardness ~ p can form a 
color field comoving with the quark pair. We plan to 
study the detailed predictions of our scenario for high p± 
quarkonium production in a future publication. 

We shall also not discuss here the special features of 
quarkonium production which a ppe ar at high xf, and 
may be related to intrinsic charm |27f| and scattering from 
light constituents p8|. Thus our discussion is limited to 
the bulk of the charmonium cross section only, which 
(at fixed target energies) originates from partons with 
(x) ~ 0.1 and (p±) ~ ra- 
il. QUALITATIVE FEATURES OF THE DATA 

The data on quarkonium production shows many inter- 
esting features and regularities. Several of them are left 
unexplained (some are even contradicted) by the dynam- 
ics assumed in the Color Octet Model (COM) and the 
Color Evaporation Model (CEM). Here we wish to make 
the phenomenological case for the three general features 
(A - C) of the production dynamics that we listed in 
Section Q. 

A. Early Color Neutralization 

Heavy quarks are produced in QQ pairs of (transverse) 
size ~ 1/m, where m is the quark mass. The pair is thus 
initially much smaller than the Bohr radius of quarko- 
nium bound states, which is of order l/(a s m). If the 
quark pair ceases to interact with its surroundings (in 
particular, turns color singlet) while it is still in such 
a compact configuration then the production rates of all 
nS states are proportional to |i?„(0)| 2 , the square of their 
wave functions at the origin. Analogous proportionality 
holds for the other 2S+1 Lj quarkonium states. 

The above argument requires no assumption about 
how the compact pair is produced. The '-R(O) propor- 
tionality test' is thus a good indication of whether the 
color neutralization occurs early or late, as measured by 
the size of the quark pair. This test should moreover be 
quite sensitive since the higher radial excitations have a 
mass near open flavor threshold. Late scattering of the cc 
system will thus affect the tp' (44 MeV below DD thresh- 
old) more than the ip (630 MeV below threshold). This 
is supported by the observation that the a(t/j')/a(J/t/j) 
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ratio is significantly reduced in central nucleus-nucleus 
(SU) collisions |^{|, as would be expected due to late in- 
teractions with comoving nuclear fragments (or plasma) . 

It has been pointed out |[o| that due to the moderate 
mass of the charm quark the wave function of charmo- 
nium is probed beyond its origin. In particular, since 
the diffractive QQ photoproduction amplitude is propor- 
tional to the square of the transverse QQ separation, the 
overlap integral between the quark pair and the bound 
state wave functions gets a negative contribution beyond 
the first node in the radial wave function of the ip' . The 
predicted ratio a(ip')/a(J/ip) = 0.17 is thus smaller 
than the |i?(0)| 2 ratio and agrees with a recent HERA 
measurement 0.150 ± 0.027 ± 0.018 ± 0.011 §|. 

The o~(ip')/a( J/ if)) ratio is remarkably universal in in- 
clastic hadroproduction processes, being nearly indepen- 
dent of the nature of the beam hadron and the target 
nucleus, and also of the energy and of xf 24 2^j3^]. This 
also holds for T(nS) states. The measured ratio for the 
directly produced J/ip cross section (from which decay 
contributions have been subtracted) moreover is consis- 
tent with |E1 



q(y/) _ r(^- 



Ml, 



0.24 



(1) 



The hadroproduction ratio ([j]) is somewhat larger than 
the one measured in diffractive photoproduction, indicat- 
ing that in hadroproduction the inelastic cross section is 
more closely proportional to the wave function at the 
origin^. 

Interesting subtleties aside, the data clearly suggests 
the relevance of the perturbative QQ wave function for 
quarkonium photo- and hadroproduction. In the COM 
and CEM approaches, on the other hand, the heavy 
quark pair turns into a color singlet only after it has ex- 
panded to a size comparable to that of the bound state. 
There is then no reason to expect the cross section ratio 
to satisfy Eq. ([!]) (although this value is also not excluded 
by those models) . We believe that the agreement of the 
quarkonium cross section ratios with expectations based 
on the wave function of the quarks created in the hard 
subprocess is not an accident. This implies that the pair 
decouples from its environment while it is still compact 
(except in the presence of nuclear comovers [fl3|). 



^In the approach discussed in this paper, the scattering am- 
plitude is proportional only to the first power of the QQ sepa- 
ration. There is also high data from the Tevatron M which 
indicates that the a(ip') /a(J/ip) ratio is still larger than in 
Eq. ([j]). This may imply an even more pointlike production 
dynamics at large p± . 



B. Reinteraction with a Color Field 

Quarkonium data provides two indications that a 
rescattering of the heavy quark pair with a comoving 
color field is important in hadroproduction. 



1. Photoproduction 

Photons do not radiate gluons^j. At the early stages of 
heavy quark creation through the 717 — * QQ subprocess 
we should therefore expect no comoving color field in the 
photon fragmentation region. With the rescattering pro- 
cess of Fig. lc thus eliminated, the production process 
should be dominated by the Color Singlet Mechanism 
(CSM) of Fig. la. It is indeed one of the remarkable facts 
of quarkonium production that the CSM works very well 
for inelastic J/ip photoproduction 33 - 35 1 , whereas the 
same model underestimates the hadroproduction cross 
section by an order of magnitude This suggests 

that hadroproduction dynamics is coupled to initial par- 
ton bremsstrahlung. 

The COM parameters which fit J/ip hadroproduction 
tend to overestimate the photoproduction cross section 
|Q,|^j35), although it is possible that the discrepancy 
could be due to higher order effects [p6|. 



2. Nuclear Target Dependence 

Cross sections of hard incoherent processes on nuclear 
targets A are expected to scale like the atomic number 
of the target, cr(A) oc A a , with a ~ 1. Modifications due 
to the A-dependence of the quark structure functions are 
minor in the presently relevant kinematic range. This is 
verified by high mass lepton pair production (the Drell- 
Yan process), for which a ~ 1.00 is observed J29 ]. 

Charm quark pairs produced in the beam fragmenta- 
tion region have large Lorentz factors and expand only af- 
ter leaving the nucleus. While compact, the uncertainty 
in the energy of the cc pairs is large and they couple 
both to open charm (DD, DA C , . . .) and to quarkonium 
{J ftp, ip', ■ ■ ■) channels. In the absence of effects due to 
partons comoving with the pair one should therefore ex- 
pect the same A-dependence for open and hidden charm, 
with a ~ 1. 

Data shows that a = 1.02 ± .03 ± .02 for D/D produc- 
tion at (x F ) = 0.031 §7|, whereas a = 0.92±0.01 for J/ip 
and ip' p9| . The deviation of a from unity for charmo- 
nium appears to be independent of Eqm, and increases 
with the charmonium momentum fraction xp- This, and 
the similar A-dependence of J/ip and ip', indicates that 



* Except via higher order resolved processes which are unim- 
portant here. 
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the nuclear suppression is not due to an expansion of the 
cc inside the target. A plausible explanation is that a fur- 
ther interaction, beyond the nucleus, between the heavy 
quark pair and comoving gluons is required for charmo- 
nium formation. 

Based on the measured A-dependence of charmonium 
production an effective 'absorption' cross section o~ a bs — 
7.3 ± 0.6 mb was obtained |38| in a Glauber framework. 
This cross section is too large for a compact cc pair 
f^ , j39|]40| and should, in the present framework, be in- 
terpreted as the joint cross section of the cc pair and the 
comoving gluon field. The size of a a bs is then reason- 
able, since the gluons are at a relatively large distance 
~ from the quark pair. An analogous interpretation 
of the nuclear suppression, albeit in a different dynamical 
picture, was earlier put forward in Ref. pi]] . 

There is no evidence for a nuclear target suppression 
of inelastic J/ij) photo- and Zepfoproduction. On the con- 
trary, there is an indication |42fl of a slight nuclear en- 
hancement at xp ~ 0.7 (a = 1.05 ± 0.03), in stark con- 
trast to the strong nuclear suppression seen at this xp in 
hadroproduction p9| . Due to the absence of comoving 
gluons and the validity of the CSM we expect a ~ 1, 
The enhancement may signal a slight antishadowing of 
the nuclear gluon distribution p3| . 

The COM and CEM assume that the process which 
turns the color octet QQ into physical quarkonium is in- 
dependent of the nature of the beam and target. Both 
models expect photo- and hadroproduction of charmo- 
nium to have the same target A-dependence, which 
should moreover equal that of open charm. 



sion for bottomonium is related to the effective distance 
~ l//i between the comoving gluon field and the heavy 
quarks, which decreases as the inverse of the quark mass. 

The T(3S I ) total cross section anomaly has been ob- 
served at both fixed target fl45],fl6"[ and collider energies 
|| , ie., for 4.2 < E CM /M r £190. Similarly, the nu- 
clear target suppression of charmonium production seems 
to be independent of the projectile type (n or p) and en- 
ergy (for 150 GeV < E LAB < 800 GeV) @,|9|. The 
discrepancy between the CSM and data on direct J/tp 
production is somewhat larger at high transverse mo- 
mentum ^> m) than for the total cross section. The 
relative contribution of P-wave decays to J/ip production 
is roughly independent of p± Q . 

The above features suggest that the anomalies ob- 
served in quarkonium production are 'leading twist' in 
the quark mass m, in the total energy Ecu and in p±, 
in the sense that the effects do not vanish as inverse pow- 
ers of any of those variables. 

In the COM, the color octet contributions which ac- 
count for direct J/ip and tp' production scale by a factor 
v 2 relative to the contributions from P-wave decays, and 
thus are relatively less important for bottomonia than for 
charmonia. Thus P-wave decay contributions dominate 
T(nS) production in the COM. In particular, T(35) pro- 
duction can only be understood by assuming [ [f6[|47jl that 
it results from the decay of an (as yet undetected) higher 
lying P-wave state. 



III. A HARD RESCATTERING SCENARIO 



C. Dependence on m, Ecu and p± 

The available data shows that the 'anomalies' of 
quarkonium production are rather insensitive to the 
quark mass m — m c , nib and to variations in kinematic 
variables such as the total energy Eqm and the quarko- 
nium transverse momentum p± . 

The measured cross section of T(35), which presum- 
ably is directly produced, exceeds the CSM prediction by 
an order of magnitude f^Q, in analogy to the original 
anomaly'. The T(IS') and T(2S) cross sections are 
more compatible with the CSM, which predicts them to 
originate almost exclusively from P-wave decays. The 
situation is thus similar to that of the charmonium sys- 
tem before the P-wave contributions were experimentally 
separated and found not to account for the bulk of the 
J cross section. It would obviously be very impor- 
tant to measure the directly produced fractions of the T 
states. 

All T(nS) states have similar nuclear target A-depen- 
dence, with a = 0.962 ± 0.014 Q. The nuclear suppres- 
sion is thus smaller than that for charmonium discussed 



In this Section we address how the features (A - C) 
of quarkonium production, mentioned in Section | and 
discussed in more detail in Section ||, can be understood 
in a QCD framework. We put forward a scenario which is 
consistent with those features, and which forms the basis 



for the explicit model studied in Section [V 



A. The Basic Heavy Quark Creation Process 

The hadroproduction of heavy quarks proceeds mainly 
via the gluon fusion subprocess^] gg — ► QQ. In our ap- 
proach (as well as in that of COM and CEM, but not of 
CSM) color neutralization occurs at a time scale which is 
large compared to the time scale 1 /m of the gluon fusion 
process. This implies that the heavy quarks are produced 
(nearly) on their mass shell. 



in Section II B 2 , but still significant compared to the 
Drell-Yan case. In our approach, the smaller suppres- 



§We neglect light quark fusion qq — > QQ, which is unimpor- 
tant for the total cross section. Similarly, higher order 'gluon 
fragmentation' diagrams Q| are irrelevant for quarkonia pro- 
duced with p± < m. 
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FIG. 3. Notation for the CM amplitude $(p). The spin 
projections Ai,2 = ±1 and A, A = ±5 all refer to the 2-axis, 
taken as the direction of fei. Only one of three contributing 
Feynman diagrams is shown. 

A basic building block of our rescattering process of 
Fig. 2 is thus the amplitude <&[ 8 1 of the gg — > QQ fusion 
process shown in Fig. 3, evaluated at leading order in the 
heavy quark momentum p in the CM frame, and with the 
quark pair in a color octet state. 

This amplitude is given in Table | in terms of the spin 
(S) and angular momentum (L) of the QQ pair (see Sec- 
tion Eqs. (§) and (§) for definitions). 

The gluon fusion amplitude $M for QQ pairs in a 
color singlet state, which is relevant for the CSM, can 
be obtained from Table Q by the substitutions f a bc — * 
0, dabcTfj — > S a b5ij/N c . It has two parts, with spin and 
angular momentum S = L = and S = L = 1, respec- 
tively. The former can directly form £0 (77) quarkonia, 
while the latter couples to 3 Po,2 (xo,a) states. The \i de- 
couples since for S = L = 1 the amplitude is symmetric 
in L z and S z , whereas the Clebsch-Gordan coefficients for 
the corresponding (J — 1, J z ) state are antisymmetric in 
L z and S z (Yang's theorem). 

The color octet amplitude $I 8 1 of Table | has contri- 
butions from S = L = 0, S = L = 1 and 5 = 0, L = 1. 
The near-threshold heavy quark cross section is domi- 
nantly S = L = 0. In the CEM flj] the S = quark 
pair is assumed to turn into physical S = 1 quarkonia 
through soft, non-perturbative gluon interactions. This 
contradicts the conservation of heavy quark spin in soft 
interactions, which is believed to be a general feature of 
QCD and follows from the non-perturbative concept of 
heavy quark symmetry pq ]. 



In the COM []16| the suppression of heavy quark spin- 
flip appears as extra powers of the velocity v ~ a s of 
the bound quarks. COM production of 3 <Si states is sup- 
pressed by v 4 in the cross section, compared to the CSM. 

B. Scenario for Perturbative Rescattering 

For the quarks to form a quarkonium bound state their 
relative momentum must be of order vm <C m. In cur- 
rent quarkonium production models the quarks are cre- 
ated directly in the gluon fusion process with a relative 
momentum of this magnitude. It is thus implicitly as- 
sumed that all later interactions are soft, commensurate 
with the bound state momentum scale. 

Here we shall consider reinteractions with momentum 
transfers of 0{fx), related to the hardness scale of the 
bremsstrahlung field in the fusion subprocess. Since v <C 
fi/m <C 1 the relative momentum of the quarks created 
in the fusion process is fairly large. The rescattering is 
hard and allows perturbative spin-flip interactions. 

It might at first appear that the rescattering physics we 
have in mind is contained in the loop correction (Fig. lb) 
to the CSM lowest order process of Fig. la. However, the 
logarithmic enhancements which favor collinear and wee 
gluon bremsstrahlung are absent in the loopQ. Conse- 
quently, the loop momentum is of order m and the spa- 
tial size of the loop is of order 1/m. The loop should 
be thought of as a vertex correction to the primary pro- 
cess, not as the rescattering envisioned in point (B) of 
Section Q. 

In order to have a rescattering which is well separated 
from the primary fusion process we must assume that 
the initial state gluon radiation gives rise to a (classical) 
color field T which is comoving with the QQ pair. The 
compact color octet QQ will dominantly interact with T 
as a color monopole (massive pointlike gluon). In these 
relatively soft interactions the internal structure of the 



"Leading logarithms in hard processes originate from tree 
diagrams. It is also straightforward to verify the absence of 
logarithms directly from the loop integral. 
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pair is preserved, in particular it remains a color octet 

The color structure and spin of the quark pair can, 
however, be changed in a harder, color dipole interaction 
with the comoving field, as depicted in Fig. 2. Here we 
consider only a single such interaction, and evaluate it 
perturbatively. The main unknown in Fig. 2 is then the 
postulated color field T^{£). Quarkonium production in 
fact offers an opportunity of detecting whether such fields 
are created in hard interactions. 

The physical picture sketched above implies that the 
classical color source T^{£) of gluons with momentum £ 
should have the following properties. 

• Since the field T originates from bremstrahlung in 
the gluon fusion subprocess it is independent of the 
beam and target. 

• Since only those components of the radiated field 
which are comoving with the quark pair are rel- 
evant, the spatial distribution of T is isotropic in 
the rest frame of the QQ pair. 

• The 3-momentum exchange £ is of order p. For 
the heavy quark propagators in Fig. 2 to be nearly 
on-shell the energy component satisfies |^°| <C \£\ 
in the quarkonium rest frame. The field T effec- 
tively acts as a time-independent color source in 
the rescattering process. 

By gauge invariance we have t^T^ = 0. Hence the field 
can be expressed in terms of its transverse (A = ±1) 
and longitudinal (A = 0) components as 

r^i) = ^2e^£)r x (l) (2) 

A 

where e x is the gluon polarization vector. As we shall 
see in Section [TV] , the quarkonium data requires that 
\T°(£)\ < \T(£)\, which is equivalent to |r A=0 (^)| < 
\T\=±i(£)\ since £„£% = and \£°\ < \£\. Hence data 
suggests that the transverse components T±\ of the color 
field dominate. 

Quark pairs that in the gluon fusion process are created 
as color singlets can evolve directly into 1 S'o and 3 Po,2 
quarkonia, as in the CSM. The influence of the color field 
r on this process is disruptive. To first order in the dipole 
interaction the pair turns into a color octet, thus losing 
its overlap with the quarkonium wave function. 

Since the direct CSM process does not require any per- 
turbative rescattering, one might expect that it will dom- 
inate the contribution of Fig. 2 for, eg., \i production. 
However, direct production of P-states is suppressed be- 
cause the quarkonium wave function vanishes at the ori- 
gin. In the rescattering process the quark pair is pro- 
duced with a comparatively high relative momentum of 
O(p) and the spatial size of its wave function increases 
before the rescattering. The expansion factor p depends 
on the time interval between the gluon fusion and rescat- 
tering processes. Our model can explain the relative rates 



°f Xc2, Xci and J ftp production provided p ~ 3. The ex- 
pansion will have a smaller effect on S'-wave cross sections 
as long as the pair stays compact compared to the size 
of the quarkonium wave function. 

In the next Section we construct a simple, specific 
model for our scenario and derive its quantitative pre- 
dictions. We show how the observed fact that the J/tp 
and ip 1 are produced unpolarized requires the field T to 
be dominantly transverse. We then have two free pa- 
rameters, the strength of the transverse field T and the 
spatial expansion parameter p, which as mentioned above 
is related to the time interval between the quark creation 
and rescattering processes. 

IV. QUARKONIUM CROSS SECTIONS AND 
POLARIZATION 

A. The Quarkonium Production Amplitude 

As illustrated in Fig. 2, the production amplitude can 
be viewed as a transition from a color octet QQ state, 
created in the gluon fusion process and described by the 
wave function $, to a quarkonium state specified by the 
wave function \F The transition is mediated by a gluon 
exchange 1Z with the color field T. The transition ampli- 
tude can be expressed as an overlap integral, 

L Z ,S Z AA, era 

x^[(p)K x - x , aS (p, q )^(q)* (3) 

which can readily be derived starting from the usual 
Feynman formulation. We work in a non-relativistic ap- 
proximation where all quark lines are on-shell and have 
3-momenta much smaller than the quark mass m. The 
z-axis is taken along the direction of the initial gluon 
momentum k\ in Fig. 2. The orbital angular momen- 
tum and spin components of the quarkonium state are 
denoted by L z and S z , whereas A, A and cr, a are the spin 
projections of the Q and Q along the z-axis before and 
after the rescattering, respectively. The relative momen- 
tum between the Q and Q is 2p before and 2q after the 
rescattering, while £ — Pf — Pi is the momentum of the 
exchanged gluon. In the quarkonium rest frame Pf = 
and £ = -Pi. 



1. The gluon fusion amplitude $ 

A standard calculation of the gluon fusion process 
gg — > QQ of Fig. 3 gives the QQ wave function, to first 
order in the quark CM momentum p, as 
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2Xd a (, c — 



ex ■ £2 P ■ e(0)*<5 



-A 



(4a) 



1Z 



- (ex pe 2 - e(2A)* + e 2 ■ p £1 • e(2A)*) } 

$ S = $ S ~> °' d «bcTij 5 ab S t] /N c ) (4b) 

As indicated, the amplitude <&M for a color singlet pair 
can be obtained from the color octet amplitude $[ 8 ' by a 
trivial substitution. The incoming gluons have colors a, b, 
momenta m(l,0j_,±l) and polarization vectors £i,£2- 
The quarks Q,Q of colors i,j have momenta (m,±p) 
and spin projections A, A along the z-axis. The spin po- 
larization vector e is defined conventionally by 



<5X^-S(-9-^A)y*«(-9,^) 



'2m 



(9) 



To first order in £/m, q/m this reduces to 



x + [51(1 - 2qf + iey^'xto*XA] } 



t 2 



e(±l) = (0,Tl,-i,0)/\/2 
e(0) = (0,0,0,l) 



(5a) 
(5b) 



(10) 



Lorentz transforming from the gg — * QQ CM to 
the quarkonium rest frame shifts any four-vector by an 
amount of 0(\£\/m). In particular, (£1 x e 2 ) z is shifted to 
(£1 x e 2 ) z [l + 0(£ 2 /m 2 )] . To the accuracy of our calcu- 
lation we can ignore the boost and use the CM amplitude 
(H) directly in Eq. (|). 

It is instructive to express the amplitudes (^) in terms 
of the spin S and orbital angular momentum L of the 
QQ pair. The amplitudes $ SSz of definite QQ spin are 
defined by 



$ 1 ' ±1 = ± $ 



($1 _i + $1 A /V2 

\ 2 ' 2 2 1 2 / ' 



$1 _1 - $_I I 

2 > 2 2 ' 2 



- 2 2 



(6a) 

(6b) 
(6c) 



The angular momentum components are defined via a 
partial wave expansion, 



Here ok are the Pauli matrices and x\ the spinors x+ — 

(1 0), xt = (o i). 

We note that the rescattering kernel 1Z consists of two 
parts. A fi = spin-conserving part which is of 0(1), 
and a [i — i part which also contains spin-flip and is of 
0(£/m, q/m). 

3. The quarkonium wave function ^ 
The quarkonium wave function \I/ in Eq. ([}]) is 

^h Sz (<?)* = ^* LLz («) 2^"(-«» °yPas. (g, -g)«(g, a-) 

(11) 

where vE'll^'Z) is the usual non-relativistic bound state 
wave function and the spin-projection operator Pss z is 
given by [fl) 



S VM2Z + 1)^- % cfy^T^ fc (7) ^ (pi , P2 ) = J_ £ (±\x^ 2 \ 2 \SS z )v(p 2 ,\ 2 )u( P x,\x) 



Ai ,A2 



where are standard spherical harmonics |50[ . The 
amplitudes $£f* are gi ven m Table Q (see also |47||). 



27ie rescattering kernel TZ 

The rescattering kernel 7?. has two terms, describing 
gluon scattering from the quark and from the antiquark, 



K x - x ,o S {P, q) = (2n) 3 S 3 (p -q+ -)K^(l, q) 



1 #2 - m fv 7° + 1 + m 



n 



2m VP2+™"™ Z 2 ^ + 
^^^ 2 -rn)nssMx 
f -75 (5 = 0) 

(s = i) 



rn 

m) (12) 



(13) 



+ (2nrS%p-q--)n^J£,q) (8) 



where 



Here we used pi — Pf/2 + q = (m, q), p 2 — Pf/2 — q = 
(m, —q) and terms of order q 2 were neglected in the last 
line of Eq. (O). A simple calculation gives for the 5 = 1 
case 

* L a f *{<!)* = n L Aq)^e(S z y ■ X l s a Xa (14) 

It is now straightforward to use Eqs. (Q), ( |l0| ) and ( |l4|) 
in Eq. (0). We contract over both pairs of indices A, A 
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and cr, a, and work to first order in £/m, q/m. The color 
structure of the field T in Eq. (|f0|) is made explicit by 



1 



(15) 



where d is the color index of the rescattering gluon and 
i',f are the colors of the quark and antiquark just be- 
fore the rescattering. The factor 1/V3 is from the color 
singlet bound state wave function. 
Finally, we use for S-wave states 



f d « ,t,* ( \ R ° 

f d 3 q , , s 
J (2^«^=° 



(16) 



where Rq is the value of the S- wave function at the origin. 
For P-wave states 

rf 3 



(17) 



where is the derivative of the P-wave function at the 
origin. 

B. The Si quarkonium cross section 

We find for the 3 Si quarkonium production amplitude 
Mi'S } . ,s. , = IT~==-^ [t\i6-^T d (£) x £ ■ e(S 2 )* 



-A2 /jz <;0 



V 67rm 3 



r^(£) 



e(Ai) • t8£ + e(A 2 



where the color factor is proportional to d a bc, 
1 



D" 



(18) 



(19) 



and Ai , A2 are the z-components of the spins of the in- 
coming gluons (A1.2 = ±1)- As can be readily inferred 
from the form of Eqs. ([I]) and ([l(]) the 1st and 2nd term 
of Eq. (|l^) correspond, respectively, to 

• the production of the QQ pair in an S = L = 
state, followed by a spin-flip interaction with a 
gluon from the color field T, and 

• the production of the pair in an S = L = 1 state, 
followed by a spin-conserving interaction with a 
gluon from the color field T°. 



The result for the amplitude squared, summed over the 
quark color indices i,j and averaged over and the 
gluon spin components Ai,A2 = ±1, is for a given spin 
component S z of the quarkonium 



xl {\T d (£)x£\ 2 + c(S z )£ 2 \T d (£)\ 2 } 



(20) 



where c(S z = 0) = 1 and c(S z — ±1) = 3. We have used 
D d D d' = S$(N 2 -4)/2N and the fact that the color field 
r is isotropically distributed. Thus, for instance, 

(T d (£) x £) x (T d {£) x £) y -» 

(r d (£) x £) z \ 2 ^ ^\T d (£) x £\ 2 

(21) 

The fixed target data p9|-pl[| shows that the polariza- 
tion of both the J/ ip and the ij/ is small and consistent 
with zero. We see from Eq. (gfj) that the condition for 
the (directly produced) 3 Si quarkonium states to be un- 
polarized is 



\T d (£)\ < T d (£) x 

In the following we shall therefore take 
T d {£) = 



(22) 



(23) 



which implies (cf. Eq. (g)) that the color field T is made 
of transversely polarized gluons (in the QQ rest frame), 
and that QQ pairs that form 3 Si quarkonia are created 
with S = L = 0. 

We thus have, for 3 Si production via gluon fusion, 



^TTTTT^^TU 160 2 3 R 2 \T d (£)x£\ 



£ 4 



(24) 



The cross section will involve an integral of this expres- 
sion over £. Hence the weighted average 



d 3 £ \T d (£) x £\ 
(2tt) 3 I 4 



(25) 



is the main free parameter of our model, which can be 
determined, eg., using the measured (direct) J/ip cross 
section. The tp' cross section then satisfies Eq. (P, as dis- 
cussed in Section ||. As we shall see below, the P-wave 
production cross sections are also proportional to T 2 , but 
in that case an additional parameter enters, which is re- 
lated to the length of the time interval between the gluon 
fusion and rescattering processes. 

Since the color field T arises through gluon radiation 
in the primary gluon fusion process it should be the same 
for all charmonium (as well as open charm) amplitudes. 
In the production of bb pairs, on the other hand, the 
gluo n radiation and thus also the color field parameter 
(p5J) will be different. 
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C. The P-wave quarkonium cross sections 

The P-wave production amplitudes are obtained sim- 
ilarly to the S- wave ones, by substituting Eqs. (Q), jlO| ) 
and ([ll]) into Eq. (|3|). However, since the P-wave func- 
tion vanishes at the origin it is now important to take into 
account the spatial expansion of the QQ pair between its 
creation in the gluon fusion process and its reinteraction 
with the color field T. 

The production amplitude (j^) can equivalently be ex- 
pressed in coordinate space as an overlap between the 
wave function of the quark pair (4>) and that of the 
quarkonium bound state (\I/) at the same spatial sepa- 
ration, 

M = J d 3 Xi d 3 x 2 d 3 x 3 $(a?i - x 2 ) r(x 3 ) **(£Ci - x 2 ) 

x [D(xi - x 3 ) + D(x 2 - £C 3 )] 
_ , Xi + x 2 



x exp 



(26) 



Here X\,x 2 are the positions of the heavy quarks at the 
reinteraction time, and D stands for the exchanged gluon 
propagator. Since \&(0) = for P-wave quarkonia the 
amplitude ( p6| ) is sensitive to the spatial extent of the 
quark wave function <f>(cci — x 2 ) at the time of reinterac- 
tion. 



According to our discussion in Section III B , the QQ 



pair is created in the gluon fusion process with a size 
~ 1/m and a relative momentum of O(p), which is large 
compared to the relative momentum q of the quarkonium 
bound state, fx ~ \£\ ^> |q|. In the (proper) time interval 
r before the pair reinteracts with T it thus expands a 
distance 



Alxl 



(27) 



The expansion time must be at least of the order of the 
spatial size of the comoving field, t > hence A|x| > 
1/m is comparable to (or even larger than) the initial 
size of the pair. We shall parametrize the expansion by 
a factor p > 1, and rescale the initial coordinate space 
gluon fusion amplitude accordingly, 



$(xi -x 2 ) 



1 



$ ( -(xi - x 2/ 
P 



(28) 



where the factor 1/p 3 preserves the normalization of the 
squared wave function. In momentum space, the rescal- 
ing d2^) implies 



(29) 



in Eq. (@). 

Intuitively it is clear that the rescaling of the argu- 
ment of $ by p will increase the overlap with the P- 
wave quarkonium wave function by the same factor (and 
hence result in a p 2 enhancement of the cross section). 



Conversely, the effect for S-waves vanishes in the limit 
where the quarkonium bound state radius 3> p/m. For- 
mally, the enhancement of the P-wave amplitude can be 
seen from Eq. (fla)), where the / a & c part (which due to 
charge conjugation invariance contributes to P-wave pro- 
duction) is linear in p, and thus gets multiplied by p in 
the rescaling (|29|) . 

We find for the spin triplet P-wave production ampli- 
tude, taking into account Eq. (|23 



M(Pj,J,)=pF —=—^8^ 
x [T d (£) xe] % J2 {LL z ;SS z \JJ z )e z {L z )*e\S z y 



where 



rpd r rpC rpd 

r — 2"' abc ij j'i' 



(30) 



(31) 



Here the QQ pair is created with S — 0, L = 1 and 
experiences a spin-flip interaction with the color field T. 
An explicit expression for the spin sum in the last factor 
of Eq. p|) may be found in Ref. f>l. Thus, 



V D7rm J *- 

r ^ [T d {£)x£] z (J = 0) 

x J -i^f [e*(J z ) x (T d (£) x £)] z (J = 1) (32) 

1 -^,(J z )[r d (f)xff (j = 2) 

The polarization tensor for a J = 2 system is given 
in Ref. g|. 

The P-wave amplitude A4( 3 Pj, J z ) depends on the 
(isotropic) color field T through the vector T d (£) x £. 
In the amplitude squared, averaged over the incoming 
gluon spins and over the momentum transfer £, the In- 
dependence thus enters through the same parameter as 
for S- waves, T 2 of Eq. (p5|). Hence the relative produc- 
tion rates and polarizations of all 3 Pj quarkonia are pre- 
dicted, and their cross sections can be compared to those 
of the 3 Si states in terms of the expansion parameter 
p. The ratios <j( 3 Pj, Jz)/&dir( 3 Si) are given in Table O, 
where (Jdir( 3 Si) is the tot al (un polarized) 3 Si cross sec- 
tion calculated in Section IV B| (cf. Eq. (pi|)). All cross 



^When Eq. (^3|) is satisfied, the rescattering kernel 1Z is lin- 
ear in the small quantities t/m, q/m (see Eq. ( p"o| ) with the 
p = term removed). For s Pj production, the only contribut- 
ing part in "1>' 8 ' of Eq. (B3) is the term oc f a bc, which is also 
linear. The amplitude (pG) thus arises from a quadratic term 
~ q-£/m 2 , however our linear approximation in £/m, q/m is 
perfectly justified. 
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sections satisfy <t( 3 P/ 2 J z ) = <t( 3 Pj,—J z ). The effective 
parameter r of Table O is defined as 



5 P 



R[/m 
Ro 



2.5 10" V 



6.510 



-V 



(cc) 
(66) 



(33) 



where the numerical values are from Ref. p5[ , with 
m c = 1.5 GeV and m fc = 4.5 GeV. Note that Table |] 
refers only to quarkonium hadroproduction through our 
rescattering process (c/. Fig. 2). The CSM mechanism 
may contribute significantly to 3 Pi production (see be- 
low). Table [n] is not relevant for quarkonium photopro- 
duction, to which our process does not contribute since 
r = (in the photon fragmentation region) due the ab- 
sence of gluon radiation from the beam photon. 





3p 


3p 


I 




3 P 2 


J z 





1 





2 


1 




r 


|r 








|r 2r 


^indir( 3 Si) 





l 

3 


1 


1 


1 3 
3 5 


TABLE II. 
polarizations. 


Relative 3 Pj 


cross 


sections and induced J ftp 


As can be 
satisfy 


seen from Table [F 


|, the total P-wave rates 


c-(Xo) : er(xi) : 


°(X2] 


= 1 


3 : 5 


(34) 



There is no data on cr(xco), but our Xci/Xc2 ratio is con- 
sistent with the value 0.6 ±0.3 measured in 7rA^ collisions 
^2j,^3|. The experimental ratio allows a CSM contribu- 
tion to Xc2 production which is about equal to the one 
given in Table |H| (c(Xci) is small in the CSM). There is 
no experimental information on the polarization of the 
3 Pj states. Nevertheless, it is interesting to note that we 
find a( 3 P 2 ,J z = ±2) = 0, contraryto the CSM where 
only this polarization is produced |24|] . 
The measured [^2| cross section ratiof^] 



0"(Xe2) 
0-dir(J/i>) 



= 5r ~ 1.8 ±.4 



(35) 



implies using Eq. ( |33j ) a value p — 2.7 . . . 3.8, where the 
lower value corresponds to the CSM contributing 50% 
of the Xc2 rate. The rather large value p ~ 3 of the 
expansion parameter is a consequence of the theoretical 
suppression of 3 Pj production due to the vanishing of 
the P-wave function at the origin (cf. Eq. (|33|)) and the 
fact that the measured a(x c i,2) nevertheless are similar 

tO <7 dir (J/tp). 



The radiative decays Xci.2 — * J/tp+"/ contribute [£2 23 
~ 40% of the total J ftp cross section. Since the %c 
states are polarized the indirectly produced J ftp's will 
in general be polarized as well. The polarization is con- 
ventionally parametrized in terms of a parameter A in 
the J /tp — > p + n~ decay angular distribution (we use the 
Gottfried- Jackson frame, where the z-axis in the J ftp rest 
frame is taken along the beam direction), 



da 



d cos 9 , 



oc 1 + A cos 2 9 



A = 



a(S z = +1) - <j(S z = 0) 
a(S z = +1) + ais, = 0) 



(36) 



As discussed above, the condition ( p3| ) implies that the 
directly produced J ftp's are unpolarized, Xdir = 0. A ra- 
diative decay contributes to the indirect 3 5*1 cross section 
according to [)52| 

O-indir^Sl, S z ) = Br( 3 Pj ^ 3 S 1 + 7 ) 



X ^\{JJ Z \1(J Z 



S z );lS z )\ 2 a( 3 Pj,J z ) 



(37) 

We show in Table || the J ftp polarization parameter 
\ndir( 3 Si) induced by the radiative decay of a 3 Pj state 
of given \J Z \. It happens that the induced J / tp polariza- 
tion is longitudinal (A < 0) for all the 3 Pj states which 
are produced by our mechanism. Using branching frac- 
tions Br(xcj -> J ftp + 7) of 27.3% and 13.5% for the xa 
and Xc2, respectively |50[, we estimate an overall polar- 
ization parameter^ \{J/tp) ~ —0.14. The measurements 
|l9| , p0| tend to prefer a value closer to zero. There is thus 
room for some Xc2 production via the CSM mechanism, 
which gives rise to fully transverse (A = 1) J/tp's. For 
o-csm{Xc2) ^ \o t ot{Xc2) the total X(J/tp) ~ -0.02. 

The tp' is only produced directly and in our mechanism 
is unpolarized for a color field T satisfying Eq. ( p3| ) . This 
is consistent with the experimental value |^l[ X(tp') — 
0.02 ±0.14. 



V. CONCLUSIONS 

Our motivation for investigating the reinteraction sce- 
nario of Fig. 2 for quarkonium hadroproduction was due 
both to regularities in the data and to shortcomings of 
alternative mechanisms, as explained in Section |l]. Here 
we shall briefly comment on some aspects of our results. 

We made several simplifying assumptions, some of 
which may need to be modified in future studies. In 
particular, 



Ref. 



gives a larger value 3.4 ± 0.9 ± 0.5 for this ratio. 



^The inclusion of the unpolarized tp' — > J ftp contribution 
does not change the numerical value of A significantly. 



11 



1. We considered only a single hard reinteraction with 
the comoving field T. 

2. We assumed T to be isotropic and independent of 
time. 

3. We assumed that only the origin of the quarko- 
nium wave function is relevant in the overlap inte- 
gral ©. 

These assumptions seem reasonable in a first attempt to 
consider the effects of reinteractions. We believe that 
further systematic studies of the environment of partons 
created in hard collisions are called for. Besides quarko- 
nium production, the flavor and azimuthal angle correla- 
tions as well as the spin dependence of open heavy flavor 
production should be informative. 

Our scenario for quarkonium production was based on 
several striking features of the data, which we in terpreted 
as the properties A - C listed in Section [A. We also 



built on the extensive experience gained from previous 
model studies of quarkonium production. Many of the 
'successes' of the present approach were thus built in 
from the start. Nevertheless, we find it non-trivial and 
interesting that so many observed features can be qual- 
itatively understood in a simple theoretical framework. 
We also found more detailed consequences of our model 
which could not be anticipated. Let us mention two suc- 
cesses and one difficulty: 

• The c(Xci)/ct(xc2) ratio is consistent with data. 
This ratio is found to be much too low, compared 
to data, in both the color singlet (CSM) (24 j and 
color octet (COM) ||(| approaches. 

• The (non-) polarization of the J/ip and Tp' . We find, 
in agreement with data, that the directly produced 
3 Si states are unpolarized (at moderate xf), pro- 
vided that the reinteraction is dominated by trans- 
verse gluon exchange, cf. Eq. (22). In both the 
CSM and COM the JAi's are produced with trans- 
verse polarization |^J,^,|2) . The color evaporation 
model (CEM) postulates that soft interactions flip 
the heavy quark spin, in violation of heavy quark 
symmetry |48[ ]. 

• Spatial expansion of the heavy quark pair. The ini- 
tially compact pair expands by a factor p before its 
reinteraction with the color field creates an over- 
lap with the quarkonium wave function. We need 
p ~ 3 to fit the observed relative rates of J/i[> and 
Xc2- Such a large expansion of the quark pair may 
be inconsistent with our approximation, which con- 
siders only the quarkonium wave function at the 
origin. 
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